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ABSTRACT 

This paper investigates the influence of magneto-centrifugally driven or simply 
magnetic winds of rapidly-rotating, strongly-magnetized T Tauri stars in causing the 
inward or outward migration of close-in giant planets. The azimuthal ram pressure 
of the magnetized wind acting on the planet tends to increase the planet's angular 
momentum and cause outward migration if the star's rotation period P* is less than the 
planet's orbital period P p . In the opposite case, P* > P p , the planet migrates inward. 
Thus, planets orbiting at distances larger (smaller) than 0.06AU(P»/5d) 2 / 3 tend to be 
pushed outward (inward) , where P* is the rotation period of the star assumed to have 
the mass of the sun. The magnetic winds are likely to occur in T Tauri stars where the 
thermal speed of the gas close to the star is small, where the star's magnetic field is 
strong, and where the star rotates rapidly. The time-scale for appreciable radial motion 
of the planet is estimated as ~ 2 — 20 Myr. A sufficiently massive close-in planet may 
cause tidal locking and once this happens the radial migration due to the magnetic 
wind ceases. The magnetic winds are expected to be important for planet migration 
for the case of a multipolar magnetic field rather than a dipole field where the wind is 
directed away from the equatorial plane and where a magnetosphcric cavity forms. The 
influence of the magnetic wind in eroding and eventually destroying the accretion disk 
is analyzed. A momentum integral is derived for the turbulent wind/disk boundary 
layer and this is used to estimate the disk erosion time-scale as ~ 1 — 10 2 Myr, with 
the lower value favored. 

Key words: stars: pre-main sequence: stars — magnetic fields: planets and satellites: 
general: accretion, accretion disks: stars: winds, outflows 



1 INTRODUCTION 

More than 250 planets have been discovered around solar- 
type stars. About 20% of them are located close to the star, 
at R < 0.1 AU, and a significant number are located as 
close as at R < 0.05 AU. According to the presently favored 
interpretation, planets form far away from the star either 
through core accretion (Mizuno 1980; Pollack et al. 1996), or 
through instabilities in the disk (Boss 2001). Subsequently, 
they migrate inward due to their gravitational interaction 
with the disk (reviewed by Papaloizou & Terquem 2006). 
For typical conditions planets migrate inward as a result of 
interaction of the planet with the disk matter. The planet 
loses part of its orbital angular momentum by overtaking 
collisions with the disk outside its orbit and it gains a smaller 
part by overtaking collisions of the disk matter inside its 
orbit. 

Some planets are expected to migrate close to the star 
where the disk properties are strongly influenced by the star 



and its rotating magnetic field. In particular, the inner re- 
gions of the disk may be dispersed as a result of heating 
by the star (Kuchner & Lecar 2002) and photo-evaporation 
(Matsuyama, Johnstone & Murray 2003). Alternatively, the 
equatorial-plane density of the accretion disk may be enor- 
mously reduced inside of what is termed the "magneto- 
spheric gap" (R < R m ) owing to to the strong magnetic 
field of the rotating protostar (Lin, Bodenheimer & Richard- 
son 1996; Romanova & Lovelace 2006). In this situation in- 
ward migration of the planet stops once the planet is some 
distance inside the gap (R < 0.63R m ) because there is no 
matter for it to transfer its angular momentum to. On the 
other hand, the magnetic interaction between a planet and 
the star may also influence a planet's migration (Papaloizou 
2007; Fleck 2008) . Indications of such interactions have been 
observed by Shkolnik et al. (2005). In these models, the 
star's magnetic field is assumed to be dipolar with sufficient 
strength that a low density cavity forms near the star. 
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In contrast, if the star's magnetic field is dominantly 
multipolar rather than dipolar (e.g., Saner 1998; Donati et 
al. 2007; Johns-Krull 2007), then there may be no cavity 
near the star. Instead the magnetic field consists of closed 
field regions and multiple open field line regions where the 
field extends to large distances. A wind can flow outward 
along the open field lines (Jardine et al. 2006; Gregory et 
al. 2006). Spectral measurements provide evidence of strong 
stellar winds from young stars (e.g., Kwan, Edwards, & Fis- 
cher 2007). 

The winds from magnetized stars may be thermally 
driven as the Solar wind is as proposed by Matt and Pu- 
dritz (2008a,b). In this case, the magnetic field still has an 
important role in the outward transport of angular momen- 
tum (Weber & Davis 1967; Matt & Pudritz 2008a,b). In 
contrast, for conditions where the thermal speed of the gas 
close to the star is small (compared with the escape veloc- 
ity) , where the star's magnetic field is strong, and where the 
star rotates rapidly, there are magneto-centrifugally driven 
winds (or more simply magnetic winds) (Michel 1969, M69 
hereafter; Belcher & MacGregor 1976. BM76 hereafter). The 
rapid rotation of the stars results of course from the accre- 
tion of high specific angular momentum from a disk. The 
magneto-centrifugal winds are driven by the star's rapidly 
rotating magnetic field rather than the thermal energy in the 
star's corona. The angular momentum and energy carried 
by the wind is extracted from the star's rotation. The mag- 
netic fields of classical T Tauri stars (CTTSs) are typically a 
thousand times larger the field of the Sun (e.g. Johns-Krull 
& Valenti 2000; Donati, et al. 2007), and these stars rotate 
much faster than the Sun with periods ~ 2 — 10 d (Bouvier 
et al. 1993). Thus magnetically driven winds may important 
for T Tauri stars. 

Magnetic stellar winds have been found in magneto- 
hydrodynamic (MHD) simulations of rapidly-rotating, disk- 
accreting stars in the "propeller" regime (Romanova et al. 
2005; Ustyugova et al. 2006). The winds are found to flow in 
opposite directions along the rotation axis transporting en- 
ergy and angular momentum (via the electromagnetic field) 
away from the star. So far simulations have been done only 
for the case of a dipole field aligned with the rotation axis. 
In this case the dipole magnetic field acts to block an equa- 
torial outflow. As mentioned for the case of multipolar fields, 
winds are expected in all directions from the star's surface. 

This paper first analyzes the torque due to magnetic 
stellar winds on close-in giant planets. This torque can cause 
the planet to migrate either inward or outward depending on 
the star's rotation period and the planet's orbital period. We 
go on to develop a model for the erosion of disks by magnetic 
winds. It is thought that close-in giant planets formed at 
larger distances (> 5 AU) and migrated rapidly inward due 
to angular momentum loss to the accretion disk (type II 
migration; Papaloizou & Terquem 2006). In addition to the 
influence of the stellar wind, a number of other effects may 
also be important in preventing the planet from migrating all 
the way into the star: (1) As mentioned above, there may be 
a magnetospheric gap near the star where the disk density 
is greatly reduced so that the inward migration stops. (2) 
The inner region of the disk may be dispered by heating by 
the star as mentioned, and it may be eroded by a magnetic 
wind. (3) There may be a strong tidal interaction which acts 
to lock the star's surface layer to the planets motion (Zahn 



1994; Marcy et al. 1997; Donati et al. 2007). Once locking 
occurs the migration due to a magnetic wind stops. 

Section 2 of this paper discusses magneto-centrifugally 
stellar winds (or simply magnetic winds) . Section 3 analyzes 
the influence of magnetic winds in pushing close-in planets 
outward or inward. Section 4 comments on the influence of 
the tidal interaction on close-in planet migration. Section 
5 investigates the role of magnetic winds in eroding and 
eventually destroying the accretion disk. Section 6 gives the 
conclusions of this work. 



2 MAGNETO-CENTRIFUGALLY DRIVEN 
WINDS 

We adopt the Weber and Davis (1976) model of an ideal 
magnetohydrodynamic (MHD) stellar wind, where the wind 
is assumed stationary and axisymmetric, and where atten- 
tion is focused on the equatorial region of the flow. In 
spherical (7?, 9, <p) inertial coordinates, the density, pres- 
sure, and flow velocity of the wind are p{R), p(R), and 
v(i?) = VR~R, + v<p(p . The magnetic field, B = BrR + B<p(j) , 
emanates from a perfectly conducting star of radius R* and 
mass M* that rotates at the angular rate f2» . It is assumed 
that p — const p 1 where 7 = const is the polytropic index. 
In the equatorial region the mass flux of the wind per unit 
solid angle is conserved, M w /4ir = R 2 pvp — const, and the 
radial magnetic flux per unit solid angle is also conserved, 
$ij/47r = R 2 B R = const. 

There are a number of conserved quantities of the flow, 
one of these gives B^/Br = (v^ — Q*R)/vr, another gives 
L — R[v4, — BrB^/^ttvr)] = const, the angular momentum 
per unit mass, and a third gives the Bernoulli constant B = 
(4 + «4)/2 + c?/( 7 - 1) - GM,/R - ^RBrB^/^pvr) = 
const, where G is the gravitational constant and c s — 
(yp/p) 1 ^ 2 i s the sound speed. The ideal MHD equations can 
be combined to give an equation for dvR/dR which can be 
integrated from the surface of the star to very large dis- 
tances for appropriately chosen parameters (BM76). Specif- 
ically, the parameters are chosen so as to allow the smooth 
integration of vr through the slow and fast magnetosonic 
singular points. 

If the star is not rotating, the radial flow velocity at very 
large distances is the Parker velocity vp = [2Cs /(7 — 1) — 
2GM,/fl,] 1/2 , where c s o is the sound speed at the initial 
radius i?* which is assumed larger than the flow speed at this 
distance. For a rotating star, one can define the reference 
velocity, 
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introduced by Michel (1969). Here, P, — 2ir/£l, is the rota- 
tion period of the star, Br, = B_r(_R«), and in the second 
line of equation (1) the star's radius is assumed to be twice 
the radius of the Sun, R* — 1.4 x 10 11 cm (Armitage & 
Clarke 1996). 

In this work we consider magneto-centrifugally driven 
winds where vm 3> vp and vm 3> c s o as discussed by M69 
and BM76. In this limit the radial wind velocity at large 
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Figure 1. Approximate radial (vr) and azimuthal (va) wind 
velocity components normalized to 2vm/3 from M69. The dashed 
line represents the azimuthal corotation velocity. 
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Figure 2. Sample magnetic field lines in the equatorial plane 
obtained using equation (4). 



distances is vm- The radial flow velocity increases from a 
small value at the star to a value (2/3)vm at the Alfven 
radius 
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These functions are shown in Figure 1. Because B^/Br = 
(w^— f2»i?) /vr, it is clear from equations (2) and (3) that this 
ratio depends only on R/Ra- Figure 2 shows the shape of 
sample magnetic field lines in the equatorial plane. For R — 




Figure 3. Sketch of a magnetic field line for conditions where the 
azimuthal wind velocity is larger than the planet's azimuthal 
velocity. 



Ra, B^/Br fa -1.55. The ratio (3 = 47rpv 2 /B 2 also depends 
only on R/Ra, and it increases with R monotonically to 
» 0.319 at R = R A . 



3 DRAG FORCE ON PLANET 

We consider a giant planet of mass M v in a circular orbit 
of radius R p ^ 1 AU. The geometry is sketched in Figure 
3. The azimuthal drag force on the planet's motion through 
the wind leads to a gradual change of the planet's angular 
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Here, T w is the timescale for the wind drag to significantly 
change the planet's orbit, A c s is the effective cross section 
of the planet discussed below, 
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is the ram pressure, and vk p = (GM, / Rp) 1 ^ 2 is the Keple- 
rian velocity of the planet. 

We consider that the effective cross-section of the planet 
A e {{ may be larger than its projected area nr p because the 
planet may have its own magnetic field or the planet may 
have an evaporatively driven wind (Vidal-Madjar, et al. 

2003) . We consider the case where the planet has a dipole 
magnetic field with surface strength B p o (Sanchez-Lavega 

2004) . Then, the effective area is 
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if this quantity is larger than irr p ; otherwise A c g — ivr p . 

The stellar wind causes outward (inward) migration of 
the planet when the azimuthal speed of the wind Vrf, is larger 
(smaller) than the planet's azimuthal speed vk p - We find 
that equality v$ = vk p holds for R p equal to a critical value, 
Rcr, that depends on the parameters of the wind, M w , Br*, 
vm, and the star's rotation rate f2* = 27r/P». For R cr <C Ra 
we find 
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Figure 4. Timescale for the wind to significantly change the 
planet's orbit T w as a function of the mass loss rate of the wind, 
M w . The bracketed quantities represent the star's surface mag- 
netic field in kG and its rotation period in days. The planet's mass 
and radius are taken to be equal to Jupiter's while the planet's 
magnetic field is taken to be 100 G. 
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which is the corotation radius of the star. A planet in a 
circular orbit of radius R cr has an orbital period of P* . This 
expression for R cr follows directly using equation (3) which 
gives U0 w Q*R (for R <C Ra) and setting this equal to the 
planet's Keplerian velocity vk p - 

Figure 4 shows the timescale T w from equation (4) 
for the magnetic wind to cause a significant change in 
the planet's orbit for different parameters. For R p > R cr 
(R p < R cr ) the wind causes the planet to move outward 
(inward). For this figure the planet's mass and radius have 
been taken to be Jupiter's and its surface magnetic field is 
taken to be 100 G. Typically, we find A e ff/7rr| ~ 2-10. The 
torque on the planet acting to move it outward is transmit- 
ted mainly by the magnetic field because we find R cr <C Ra 
where (3 = 47rpv 2 /B 2 < 1. The increase in the planet's 
angular momentum equals a corresponding decrease in the 
star's angular momentum. Assuming R cr <C Ra and using 
equation (7) we find the estimate 
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where it is assumed that A c s ^ Tvr p . In the second line the 
planet's mass and radius is assumed equal to Jupiter's and 
its surface magnetic field as 100G. The mass of the star is 
Mq and its radius is 2Rq. 

Alternatively, the torque on the planet can be thought 
of in terms of a DC current circuit following the model of Io's 
interaction with Jupiter's magnetosphere by Goldreich and 
Lynden-Bell (1969). This is applicable if the Alfven wave 
transit time from the planet to the star and back is less 
than the time for the planet to move a distance of the order 
of its diameter relative to the wind (Drell, Foley, & Ru- 
derman 1965). The torque is T = —R p IgBR(R p )(2r e g)/c, 
where Ig is the poloidal current flowing in the 9— direction 
through planet's magnetopause. This current then flows 



down the wind field line intersecting the top of the planet 
at = Tes/Rp', it then flows through the star's atmo- 
sphere; and finally it flows outward on the wind field line 
which intersects the bottom of the planet at 9 — —r c s/R P - 
The voltage drop across the planet in its reference frame is 
AV = — (2r e g )(v# — vk p )Br(R p )/c, which acts as a battery. 
For representative values of R p — 10 12 cm, Br = 20 G, 
r c ff = 7.15 x 10 9 cm (Jupiter's radius), and v$ — VKp = 10 6 
cm/s, one finds AV ~ 3 x 10 9 V. The parts of the circuit 
along the field lines have negligible resistance. Thus the to- 
tal current Ig — AV/1Z, where 1Z is the resistance of the 
parts of the circuit in the star's turbulent atmosphere and 
in the magnetopause of the planet. Setting the magnetic 
field torque (equation 4) from the MHD model equal to the 
circuit model torque T gives the total poloidal current flow, 
Ig ~ 3.5 x 10 10 A for the mentioned conditions. This implies 
that the circuit resistance is 1Z ~ 0.09 Ohm. The circuit 
picture is helpful for understanding the interaction but a 
consistent treatment requires a treatment using the MHD 
equations (Parker 1969). 

Observations and theory relevant to the magnetic cou- 
pling of hot Jupiters to their stars is discussed by Cuntz, 
Saar, & Musielak (2000), Shkolnik et al. (2005), Mclvor, 
Jardine, & Holzwarth (2006), and Preusse et al. (2006). 

The stellar wind acts to spin-down the star giving a 
torque on the star — (R A sm 2 (9))Q,AI w (Weber & Davis 
1967) as discussed by Matt and Pudritz (2008a, b) for the 
case of a thermally driven wind. Here, (R\ sin 2 (6')) 1 ^ 2 is the 
mass- weighted, angle averaged Alfven radius for the angular 
momentum outflow in the wind. This average Alfven radius 
may be significantly smaller than the the values Ra consid- 
ered above. The angular momentum of the star is feM,i? 2 Sl, 
with k « 0.2 for stars younger than 10 7 yr (Armitage & 
Clarke 1996) so that the spin-down time-scale of the star is 



0.2M,Ri 
M W (R 2 A sin 2 



(9) 



For the reference parameters of equation (8) and M w ~ 
10 -9 — lO _lo M0/yr, the spin-down time T ac i is less than the 
wind time T w . Under such conditions the planet would at 
first have R p > R cr and be pushed outward by the wind, 
and subsequently it would have R p < R cr and be pushed 
inward. 



4 INFLUENCE OF TIDAL INTERACTION 

We consider that initially the planet is at a relatively large 
distance with R p > R cr = (GM t /fll) 1/3 but with R p < 
Ra- We assume that the inner region of the disk (< 1 AU) 
has been expelled so that there is no migration due to the 
planet's interaction with the disk. 

For R p > R cr , the rotating magnetic wind tends push 
the planet outward with a time-scale T w given by equation 
(8). At the same time, the magnetic wind causes the star to 
spin-down with a time-scale T ac i given by equation (9). If the 
spin-down time is longer than T w the planet will be gradu- 
ally pushed outward. This case is not considered further. In 
the opposite case where T s d < T w , the planet and star will 
evolve to have R p < R cr , and the wind will tend to push 
the planet inward, again with the time-scale of equation (8). 
However, the spin-down of the star acts to increase T w and 
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very likely T w /T sd with the result that the inward motion 
of the planet will cease. 

As the planet moves inward it raises an increasing tide 
on the star. We assume that the planet itself is tidally locked 
and in a circular orbit (Marcy et al. 1997). The time-scale for 
the establishment of corotation of the entire star is estimated 
by Zahn (1994), Marcy et al. (1997), and Donati et al. (2008) 
as 

t ~~ "^UfcHife)'- (io) 

where we have assumed a solar mass star of radius 2Rq. 
The values of T tc ii are probably longer than the age of the 
systems except for very massive planets. 

Marcy et al. (1997) and Donati et al. (2008) suggest that 
only the convective envelope of the star of mass 5M* <C M, 
is brought into corotation or "locked" to the planet. The 
time-scale for this to occur is reduced from Ttdi by a fac- 
tor ~ (<5M*/M*) 2 . The core of the star then rotates slower 
than its surface. For a Jupiter mass planet, R p = 10 12 cm, 
and <5M*/M* < 0.005, we have T t( u < T w . Of course, once 
the star's surface becomes locked to the planet, the effect of 
the magnetic wind on the planet's orbit vanishes. With the 
planet /stellar surface locked, the loss of angular momentum 
from the planet/steller surface system due to a magnetic 
wind acts to reduce the planet's period if it is initially rela- 
tively large. However, the reduction of the angular momen- 
tum is limited for a fixed <5M*. This limit may such that 
a sufficiently massive planet will fall into the star. For a 
less massive planet the system behavior for further angular 
momentum loss is unclear. 

Note that the influence of the wind on the planet's 
orbit decreases as M p increases and as the star's rota- 
tion period increases (equation 8). On the other hand the 
tidal interaction increases as M p increases, and it increases 
strongly as R p decreases (equation 10). The remarkable ob- 
ject r— Bootis has a planet mass M p » 7.5Mj and period 
P p — 3.3 d and is locked to the rotation of the star's surface 
(Donati et al. 2008; Catala et al. 2007; Butler et al. 1997). 
On the other hand the famous object 51— Pegasi (Mayor & 
Queloz 1995) has a planet of mass w 0.7Mj and a period 
of 4.23 d which orbits 51— Peg which has a rotation period 
P, « 37 d. 

Consider what may happen to a giant planet after its 
fast inward migration ends at a distance say R p i < 0.1 AU 
owing to the planet entering a magnetospheric gap or the re- 
gion of the disk which has been blown away by a stellar wind. 
The planet's distance R p may be unaffected by the star's 
wind owing to the star's rapid spin down. If the planet's 
mass is large enough relative to the mass of the convective 
envelope, then it can cause tidal locking of the star's surface 
as observed in r— Boo. For a smaller mass planet there is no 
locking as is the case of 51— Peg. Another possibility is that 
R p changes due to the influence of the magnetic wind. The 
spin-down of the star is expected to give P p < P* in which 
case the wind causes a slow inward motion of the planet or 
decrease of R p . The planet moves inward until there is tidal 
locking with the star's surface (the time-scale for which is 
oc R p ) The tidal locking eliminates the wind affect on the 
planet's orbit. The migration due to the stellar wind allows 
smaller mass giant planets to produce tidal locking. 



For tidal locking in a given period T tc u, the distance 
where this occurs is R p (tdl) oc R, (M P /5M*) 1//3 . For stars of 
decreasing mass M» < M© (but M* > 0.1M©), the star's 
radius decreases strongly, and the mass of the convective 
envelope <5M* increases rapidly (the star is fully convective 
for M* < (0.3 - 0.4)M Q ) (Chabrier & Baraffe 2000). Thus 
R p (tdl) decreases rapidly as M v decreases from M©. 



5 EROSION OF THE DISK BY THE STELLAR 
WIND 

Stellar winds have long been considered as a mechanism for 
the dispersal of accretion disks (Cameron 1973). Theoret- 
ical and simulation studies of this process have been out 
assuming an entirely hydrodynamic interaction (e.g., Canto 

6 Raga 1991; Richling & Yorke 1997; Hollenbach, Yorke, & 
Richstone 2000; Soker 2005). Here we discuss the interaction 
of a high velocity magnetized stellar wind with the accretion 
disk. The considered geometry is shown in Figure 5. The 
presence of the magnetic field leads to Reynolds numbers 
which are sufficiently large that the wind/disk boundary 
layer is strongly turbulent. In the hydrodynamic case the 
Reynolds numbers are much smaller so that the boundary 
layer is probably laminar (see Schlicting 1968, ch. 7). In the 
following we derive an estimate of the mass loss rate from 
the disk. 

It is useful to consider the physical conditions in the 
wind at a distance R = 1 AU for fiducial conditions of wind 
density n w — 10 5 cm~ 3 , wind speed v w = 500 km/s, time av- 
eraged magnetic field B w =0.1 G (predominantly toroidal), 
and ion (proton) and electron temperatures Ti = T e = 10 5 
K. At this and larger distances note that (i) the wind ve- 
locity is predominantly radial, (ii) it is super fast magne- 
tosonic, and (iii) it is much larger than the Keplerian veloc- 
ity of the disk matter. The ion and electron gyro- frequencies 
are uj ci « 960 s _1 and cj ce w 1.8 x 10 6 s _1 ; the ion and 
electron gyro-radii are r gi k 3 x 10 3 cm and r ge w 70 
cm; the ion and electron collision times (oc T 3,/2 /n) are 
Ti w 470 s and r e ~ 11 s; and the ion and electron mean- 
free paths are ii = l e « 1.4 x 10 9 cm (Braginskii 1965). 
Thus we have LU ci Ti « 4.6 x 10 5 and ui ce r e R! 2 x 10 7 . In 
the absence of a magnetic field, the kinematic viscosities 
of ions and electrons are u oi « v 2 hi Ti ^4x 10 15 cm 2 s _1 
and foe ~ v 2 he r e ~ 1-7 x 10 17 cm 2 s _1 , where v t hi,e are 
the ion or electron thermal speeds. Thus without a mag- 
netic field the Reynolds number, Re — Rv w /foe ~ 4400, 
is such that a boundary layer flow is laminar. With the 
magnetic field included, there are five different viscosity co- 
efficients. However, for the considered problem the impor- 
tant viscosity coefficient is that for momentum transport 
across the magnetic field; for example the momentum flux- 
density component YIb.9 = —pv±R~ 1 dvR/d9. For the ions 
it is w±i ~ voij '{uJciTi) 2 ~ 1.9 x 10 4 cm 2 s _1 and for elec- 
trons z/± e « vo e /(Lo ce T e ) 2 « 440 cm 2 s _1 . Using the vis- 
cosity v±i, the effective Reynolds number for the wind is 
Re w = Rv w /v_n ~ 5 x 10 17 so that the boundary layer 
flow is strongly turbulent. This is in contrast with the non- 
magnetic case where the boundary layer flow would be lami- 
nar. The large reduction of the viscosity results from the par- 
ticle step size between collisions being a gyro-radius rather 
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Figure 5. Sketch of disk- wind turbulent boundary layer with the 
outer oblique shock, n is the normal to the shock wave, v w is the 
wind velocity, h is the disk half-thickness, and H is the vertical 
height of the shock. 



than a mean-free path. Thus the estimated Reynolds num- 
ber also holds for a turbulent magnetic field. 

Similarly, the important heat conductivity coefficient 
is that for the heat flux across the magnetic field, say qg — 
-K±R~ 1 d(k B T)/dO, where k± ta K ±i w 2nv 2 hi n / (uj ci Ti) 2 w 
3.8 x 10 9 (cm s) _1 and where ks is Boltzman's constant. For 
this heat conductivity, the heat flow from the wind into the 
disk is negligible compared with the energy flux per unit area 
from the disk 3GM t M a /(8nR 3 ) for accretion rates M a = 
KT 10 - lO^Moyr" 1 . 

At the place where the wind encounters the much denser 
disk, there is a weak oblique shock at a height H(R) <C 
R above the equatorial plane and above the disk which 
has a half-thickness h(R) < H(R) as sketched in Figure 
5. The angle between the incident flow and the shock is 
/3 « Rd(H/R)/dR <C 1. In general we expect to have /3 > 0. 
In passing through the shock the flow is deflected through 
an angle S — 2/9/(1 + 7) = 3/3/4 (for 7 = 5/3) away from 
the equatorial plane, and the flow speed is reduced by a 
small fractional amount. The region between h and H is re- 
ferred to as the boundary layer. The influx of wind matter 
into this layer is — dS ■ (p w v w ) ~ dSp w v w R[d(H /R)/dR], 
where dS ~ RdRd(j> is the area element of the shock on the 
top side of the disk. The Keplerian velocity of the disk is 
small compared with the wind velocity for R ^ 1 AU, and 
it is neglected. The density in the boundary layer varies 
from p(R, h) 3> p w at the surface of the disk to p w at 
z — H. The time-averaged radial flow velocity varies from 
\vr(R,K)\ <C v w to vr = v w at z = H. Here, we neglect 
v R (R, h). 

For stationary conditions, the conservation of mass and 
radial momentum in the annular region abed of the boundary 
layer on the top side of the disk - sketched in Figure 6 - gives 



^(RFr) 

Ir^ 



= R 



= R 



R 



d(H/R) 



R 



dR 
d(H/R) 



dR 



Pw "t - 



PwV u 



1 dM d 
2tt dR 



(11) 
(12) 



where the term on the left-hand side is from the vertical 
sides of the region and the right-hand side is due to the 
sides ab and cd. Here, Fr is the mass flux and Fr is the 
radial momentum flux both per unit circumference of the 
top side of the disk. That is, 



II 



dz (pv R ) 



F p — 

R — 



r 

Jh 



dz (pv%) 



(13) 



where the averages are over the turbulent fluctuations. The 
mass loss rate of disk per unit radius due to entrainment 
is dMd/dR. The disk matter influx to the boundary layer 
brings in negligible radial momentum. 

We now integrate equations (11) and (12) from an inner 
radius _Ri n , where F^ and Fr are negligible, to distance R. 
We define an average radial velocity in the boundary layer as 
u(R) = J dz(pv R }/ J dz(pv R ). Thus we obtain an equation 
for the mass loss rate from the top and bottom sides of the 
disk between R — R[ n and R, 



M d = 4tt 



RdR 



R 



d(H/R) 



dR 



Pw^u 



(14) 



which is von Karman momentum integral for this prob- 
lem (Schlicting 1968, ch. 8). Using the fact that p m v w = 
we can evaluate this integral at the outer ra- 
dius of the disk Rout- This gives the mass loss rate of the 
disk, 



(15) 



The main unknown in this equation is u. 

The average velocity u depends on the vertical profiles 
of density and radial velocity which are not known. We ex- 
pect the profiles, for example (vr(z)} to be substantially 
different from those of laboratory turbulent boundary lay- 
ers over solid surfaces (e.g., Schlicting 1968, ch. 23; Roy & 
Blottner 2006). The main reason for the difference is that 
the density at the surface of the disk p(R, h) is many or- 
ders of magnitude larger than the wind density p w . For a 
laboratory boundary layer, a mixing-length model of the 
momentum transport gives (z') 2 (d(vR) / dz') 2 = const (with 
z = z — h) , and this gives the well-known logarithmic veloc- 
ity profile (see e.g. Schlicting 1968). For this profile most of 
the change of velocity is quite close to the wall (z' — 0). In 
contrast, for the disk boundary layer a mixing length model 
gives p(z')(z') 2 (d{vR) /dz') 2 = const. Because of the density 
dependence, the change in the velocity occurs relatively far 
from the wall. An important consequence of this is that the 
average velocity u is much smaller than v w . If we make the 
weak assumption that u is less than v m /2, then we have a 
lower limit on the disk mass loss rate Md J> (H/R) out M m . 
For an initial disk mass of M d = 0.02Mq (e.g., Kuchner 
2004), (H/R) out = 0.2, and M w = 10" 9 M Q yr _1 , the disk 
loss time is T d <; M d /M d = 10 8 yr. 

However, the abovementioned discussion suggests that 
u is much less than v w . Thus the above upper limit on T d is a 
large over estimate of the erosion time. Because u is necessar- 
ily larger than the local escape velocity v csc = (2GM» / R) 1 ^ 2 
in order for the matter flow in the boundary layer to become 
unbound from the star. Thus we find an upper bound on the 
mass loss rate from the disk due to the wind 



M A 
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< 17 



500km/s J V50AU 



Rout \ / (H/R)out 



V 0.2 



M w ,(16) 



assuming v esc <C v v 
0.02M Q , {H/R) out : 



For an initial disk mass of Md 
0.2, Rout = 50 AU, and M w 
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h Hz 

Figure 6. The right-hand part of the figure shows the control 
volume abed discussed in the text. The left-hand side of the figure 
shows the vertical profiles of the density p and radial velocity vr 
in the boundary layer. Here, pd is the density at the surface of 
the disk, and at z = H, vr = v w . 



Compared with the magneto-spheric cavity model pro- 
posed for halting the inward migration of planets (Lin et al. 
1996; Romanova & Lovelace 2006), the influence of the mag- 
netic wind may allow for a broader distribution of distances 
of exosolar planets (Pont 2007). 

The influence of the magnetic wind in eroding and even- 
tually destroying the accretion disk is analyzed. A momen- 
tum integral for the boundary layer is derived and used to 
estimate the disk erosion time-scale Td- For the considered 
conditions we find Td ~ 1 — 10 2 Myr with the lower value 
favored. 

We thank Ralph Pudritz for valuable criticism of an 
earlier version of this work. This work was supported in 
part by NASA grants NAG5-13220 and NAG5-13060 and 
by NSF grant AST-0507760. 



10" 9 M Q yr" 1 , the disk loss time is T d > M d /M d = 1-2 x 10 6 
yr. Equation (16) has a different dependence on H/R and 
gives different values compared with the hydrodynamic es- 
timates (Hollenbach et al. 2000). 



6 CONCLUSIONS 

We conclude that magneto-centrifugally driven winds of 
strongly magnetized, rapidly rotating T Tauri stars may af- 
fect the orbits of close-in giant planets that are inside the 
Alfven radius of the wind. The magnetic wind may of course 
disperse the inner part of the disk and thereby halt the in- 
ward migration of a giant planet. However, this work focuses 
on the affect of the azimuthal ram pressure of the magne- 
tized wind on the planet. This gives a torque which tends to 
increase (decrease) the planet's angular momentum thereby 
causing it to move outward (inward) if the planet is at a 
distance R larger (smaller) than the corotation radius of the 
star, R cr = (GM./nlfl 3 ~ 0.07AU(P»/5d) 2 '' 3 for a solar 
mass star, where Q* = 2n/P* is the star's angular rotation 
rate. Such winds are likely to occur in T Tauri stars where 
the thermal speed of the gas close to the star is small com- 
pared with the escape velocity, where the star's magnetic 
field is strong, and where the star rotates rapidly. The mag- 
netic winds are expected to be important for giant planet 
migration for cases where the star's magnetic field is a mul- 
tipolar field rather dipolar. For an approximately aligned 
dipolar magnetic field the wind is directed away from the 
equatorial plane and a magnetospheric cavity forms. (Lin 
et al. 1996; Romanova & Lovelace 2006). We find that the 
time-scale for the magnetic wind to change the planet's orbit 
T w ranges from 2 — 20 Myr for stellar rotation periods of 3 — 5 
d and surface magnetic fields of 1 — 3 kG. This time-scale is 
of the order of or longer than the estimated spin down time 
of the star T sd due to the magnetic wind. 

The migration of a close in giant planet due to the star's 
wind may be strongly affected by the tidal interaction. The 
time-scale for the planet to establish tidal locking T t di of the 
entire star is estimated to be longer than the age of most 
systems. However, locking of the convective envelope of the 
star may occur on a much shorter time-scale. This time scale 
varies as the sixth power of the distance to the planet. Once 
locking of the stars surface occurs, the influence of the wind 
on the planet's orbit vanishes. 
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